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Abstract 



Various features of domain walls in supersymmetric gluodynamics are dis- 
^ . cussed. We give a simple field-theoretic interpretation of the phenomenon of 

^ I strings ending on the walls recently conjectured by Witten. An explanation 

of this phenomenon in the framework of gauge field theory is outlined. The 
phenomenon is argued to be particularly natural in supersymmetric theories 
which support degenerate vacuum states with distinct physical properties. 
The issue of existence (or non-existence) of the BPS saturated walls in the 
theories with glued (super)potentials is addressed. The amended Veneziano- 
Yankielowicz effective Lagrangian belongs to this class. The physical origin 
of the cusp structure of the effective Lagrangian is revealed, and the limita- 
tion it imposes on the calculability of the wall tension is explained. Related 
problems are considered. In particular, it is shown that the so called discrete 
anomaly matching, when properly implemented, does not rule out the chirally 
symmetric phase of supersymmetric gluodynamics, contrary to recent claims. 



1 Introduction 



Recently there have been a renewed interest in the study of theoretical aspects 
of = 1 supersymmetric (SUSY) gauge theories. In addition to the calculation of 
exact effective potentials and conjectured dualities between theories with distinct 
gauge groups [3, it has been realized that in some supersymmetric theories there 
exists a class of dynamical objects whose energy can be calculated exactly 0. Those 
are the domain walls interpolating between discrete vacua which are typical for many 
SUSY gauge theories. The remarkable fact is that the energy density (tension) of 
these walls is exactly calculable even in the strong coupling regime. 

For supersymmetric gluodynamics, the theory of gluons and gluinos with no 
matter, the calculation of the energy density was carried out in Ref. in an 
indirect way. The key ingredient is the central extension of the = 1 superalgebra, 

MQ^} = ^ (^)„, / rf^xV(TrA% (1) 



where Q]^ is the supercharge, A is the gluino field, and {a)^^ = {a^, —i, — cr^}^^ 



IS 



a set of matrices converting the vectorial index of the representation (1,0) of the 
Lorentz group in the spinorial indices. The commutator ([^) is given for SU{N) 
gauge group; the parameter N stands for the number of colors. The right-hand side 
of Eq. (|1|) is a reflection of "geometric" anomalies of SUSY gluodynamics (i.e. that 
in the trace of the energy- momentum tensor plus its supergeneralizations) . 

The integral over the full derivative on the right-hand side is zero for all localized 
field configurations. It does not vanish, however, for the domain walls. Equation 
(0) implies that the tension of the domain wall is 

^ (TrA2)oo-(TrA2)_^ , (2) 



^7r2 



where the subscript ±oo marks the values of the gluino condensate at spatial in- 
finities (say, aX z ±oo assuming that the domain wall lies in the xy plane). The 
existence of the exact relation (0) is a consequence of the fact that the domain 
wall in the case at hand is a BPS-saturated configuration preserving 1/2 of the 
original supersymmetry. A general discussion of BPS saturated domain walls was 
given in 0], while a particular wall realization in the framework of the amended |^ 
Veneziano-Yankielowicz effective Lagrangians [^] were studied in some detail in ^ 
and §. 

On the other hand, a theory related to supersymmetric gluodynamics was an- 
alyzed recently from the point of view of D-brane physics 0. In this picture the 
domain walls also appear naturally. Moreover they seem to have some rather sur- 
prising properties. These properties are natural from the D-brane perspective but 
were considered unusual (even paradoxical) from the field-theoretical point of view 
ip. One of such features is an "abnormal" A^ dependence of the wall tension. The 
wall energy density of some BPS saturated walls scales as A^, rather than A^^^ ^ 
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dependence one might expect apriori from glueball solitons. The second surprise 
is that the confining QCD string emanating from the probe color charges (quarks) 
on one side of the wall can terminate on the wall, without penetrating on the other 
side. 

So far these features had no satisfactory explanation in the field-theoretical 
framework. One of our tasks is to understand how this works in field theory, at 
least at a qualitative level. We show that both aspects - the N dependence and 
termination of the flux tubes on the walls - are quite natural consequences of the 
pecuhar gauge dynamics. We suggest various toy models, which are simpler than 
SUSY gluodynamics but still carry essential features of the phenomena under discus- 
sion, to substantiate our qualitative observations. For instance, an Abelian model 
is presented where probe fractional charges generate induced "mirror" fractional 
charges on the wall in spite of the fact that at the fundamental level the model 
contains no fractionally charged fields. 

Then we turn to the Veneziano-Yankielowicz effective Lagrangians. Previously 
they were exploited as a framework for quantitative analysis of the domain walls 
1^. A BPS wall interpolating between a chirally symmetric vacuum |^ and one 
of the conventional vacua with (AA) 7^ was explicitly built. However, building a 
wall interpolating between two neighboring chirally asymmetric vacua turned out 
to be a much harder endeavor. The task had not been solved in Ref. 0. Moreover, 
it was argued later ||^ that such walls do not exist within the framework of the 
Veneziano-Yankielowicz effective Lagrangians. 

A crucial feature of such a Lagrangian emerging in SUSY gluodynamics is dis- 
continuity of the superpotential W. The realization compatible with all symmetries 
of the underlying theory requires, with necessity, a "glued" potential P, with 
distinct sectors and matching lines along the boundaries of the sectors. We explain 
the physical nature of this phenomenon. The sector pattern, with cusps, reflects a 
restructuring of heavy degrees of freedom (which were integrated out) in the process 
of an adiabatic variation of the light degrees of freedom. A level crossing takes place 
in the heavy sector of the theory. Precisely for this reason, one cannot construct the 
domain wall from the effective Lagrangian if the wall crosses the cusp. The presence 
of the cusps prevents one from being able to use this potential for calculating a wall 
proflle if the fleld conflguration along the wall crosses the cusp somewhere in space. 
In fact, if one naively tries to do this in the presence of the cusp an apparent paradox 
arises - the wall in the effective theory seems to have a lower energy density than 
the BPS bound on this quantity in the original theory. The missing energy density 
is contributed by the excitation of the heavy modes which are necessarily excited 
when the light flelds take values in the vicinity of the cusp. 

The statement is thoroughly illustrated by two toy models. The phenomenon is 
quite general and may be considered in supersymmetric as well as non-supersymmetric 
context. 

The chirally symmetric vacuum, (AA) = 0, is inherent to the Veneziano-Yankielowicz 
Lagrangian. Recently it was claimed [^] that a discrete anomaly matching rules 
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out the existence of such phase, at least its most straightforward reahzation. We 
make a digression to show that the claim is due to an inconsistent treatment of the 
discrete anomaly matching. In SUSY gluodynamics and similar theories the discrete 
anomaly matching imposes no constraints on the spectrum. The only information 
one gets are rather mild constrains on certain amplitudes following from the clas- 



sical symmetries that become anomalous at the quantum level (see e.g. |TT1)- The 
existence of a discrete anomaly-free subgroup adds no new information. 

The organization of the paper is as follows. Section 2 is devoted to general 
aspects of the BPS walls in SUSY gluodynamics. We analyze, qualitatively, the N 
dependence of the wall tension and visualize the phenomenon of strings ending on 
the domain walls. An analogy between the walls in SUSY gluodynamics and the 
axion domain walls in gauge theories with monopoles is worked out. 

In Sect. 3 we turn to a more quantitative discussion based on the amended 
Veneziano-Yankielowicz effective action. The issue of glued potentials is studied 
here. In Sect. 3.3 we consider an explicit example of a supersymmetric theory which 
upon integrating out heavy modes generates an effective potential with cusps. In 
this model we calculate explicitly the cusp contribution to the wall tension and show 
how the apparent contradiction with the BPS bound is resolved. In Sect. 3.4 aspects 
of the general theory of glued (super)potentials are presented. 

In Sect. 4 we discuss other domain walls in gauge theories obtained as a Kaluza- 
Klein reduction on topologically non-trivial space-time manifolds. The specific ex- 
ample considered refers to i?4 x 5*1. An unconventional dependence of the wall 
tension arises which may be related to D-branes. 

Section 5 is devoted to the issue of the discrete anomaly matching and the chirally 
symmetric vacuum of SUSY gluodynamics. 

Finally, Sect. 6 contains a brief summary and discussion of our results. 



2 N Dependence, Flux Tubes Ending on the Walls 
and All That 



2.1 SUSY Gluodynamics 

We consider the supersymmetric generalization of pure gluodynamics - i.e. the 
theory of gluons and gluinos. At the fundamental level the Lagrangian of the model 



has the form [|T2| 
£ = 



1 ^ ~ 1 



JO 



(3) 



where the spinorial notation is used. In the superfield language the Lagrangian can 
be written as 



C = ArTr / d^eW^ + H.c. , 
4^2 J 



(4) 
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where 

1 _ 1 id 

Here d denotes the vacuum angle. Our conventions regarding the superfield for- 



mahsm are summarized e.g. in the recent review |T^. We will limit ourselves to 
the SU (N) gauge group (the generators of the group T'^ are in the fundamental 
representation, so that Tr (T'^T'') = (1/2)5'"''). 

SU{N) supersymmetric gluodynamics has a discrete symmetry, Z2N, a (non- 
anomalous) remnant of the anomalous axial symmetry generated by the phase rota- 
tions of the gluino field. The gluino condensate (AA) is the order parameter of this 
symmetry. The discrete chiral symmetry may or may not be spontaneously broken 
0. Therefore, there exists a set of distinct vacua labeled by the value of the gluino 
condensate. In the phase with the broken chiral symmetry Tr(AA) = exp{2TTik/N) 
where = 0, 1, 2, — 1 (for = 0). In the chirally sjTiimetric phase Tr(AA) = 0. 
The field configurations interpolating between different values of (AA) at spatial in- 
finities are topologically stable domain walls. Although the theory is in the strong 
coupling regime one can derive an exact lower bound on the surface energy density 
(tension) for such a wall ^ 



e > 



N 
8^ 



(TrA^)oo-(TrA^)_oo • (5) 



In our normalization the condensate (Tr A^) scales as A^ in the large A^ limit. 

One may consider two types of walls. The wall of type I connects a vacuum with 
the spontaneously broken chiral symmetry with the symmetric vacuum ((TrA^) = 
0). For such a wall the BPS bound for the tension is 

e ~ 0{N'^) . 

The walls of type II connect two adjacent (or close) chirally asymmetric vacua, e.g. 
k = and k = 1 (or k = and k = 2, etc.). Even though for each of these vacua 
the order parameter is of order A^, the difference between the order parameters is 
0(1). The BPS bound for the tension, therefore, is 

e ~ 0{N) . 

Let us assume for the moment that the BPS-saturated walls of type II do indeed 
exist in SUSY gluodynamics. Although we cannot prove this at present, there are 
no visible reasons forbidding them [|. 

The question then arises as to how one can understand the large A^ scaling of 
the wall energy density from the point of view of the effective field theory which 
describes dynamics of the low lying physical states, mesons and glueballs and their 
superpartners. 



^ Even if for some reasons we do not understand at present the type II walls would turn out to 
be not BPS-saturated, it is natural to expect that their tension is of order of the BPS bound (||). 
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At large the mesons and the glueballs should have masses of order 1, trilinear 



couplings of order 1/iV, and so on |jTj]. This is conveniently encoded in an effective 
Lagrangian of the form 

C = N^F[Mi, dMi] , (6) 

where {Mj} is a set of fields representing all relevant degrees of freedom, mesons 
and glueballs. The value of the functional F itself and all its derivatives at the 
minima should be independent of A^. This would ensure the proper A^ dependence 
of the masses and coupling constants. Now, suppose we have a solution of classical 
equations of motion Mwaii which describes a wall configuration interpolating between 
two distinct minima. Since A^^ is an overall factor in Eq. (^, at first sight one may 
expect that F[M„aii] = 0(1), and, therefore, the wall tension e = 0{N'^). Such 
a situation is standard in the soliton physics. This is perfectly in line with what 
we get for the type I walls, but is in apparent contradiction with the BPS energy 
density of the type II wall. 

Can one avoid this conclusion? The answer is yes. Consider a function F which 
is nonanalytic so that although all its derivatives are 0(1) at the minimum M^,, at 
some finite distance M = + 6M (with 6M = 0{1/N)) the derivatives become 
large (O(A^), for example). Then, if the wall solution Mwaii passes through this 
region of the field space, the standard counting leading to e = O^N"^) does not 
work. An extreme example of such a situation arises if the effective Lagrangian has a 
structure of A^ distinct sectors in the space of fields, and no single analytical function 
F exist. The A^ sector structure then is an implicit source of A^ dependence. Such a 
potential is singular (or rather has a singular first derivative) along the boundaries 
of the sectors, being glued from different pieces along the boundaries. As we will 
discuss later, this situation can arise due to the fact that the state that was the 
ground state in one sector, becomes an excited state in another, and vice versa. At 
the boundary there are degenerate states, and the level crossing occurs. Due to the 
cusps at the boundary the naive estimate of the tension presented above does not 
work in the case of the glued potential. As we will see, the effective potential in 
SUSY gluodynamics has precisely this A^-sector structure. 

Of course, in the full theory everything is smooth. One can ensure the smoothness 
of an effective Lagrangian by including more fields in it. Those extra fields will not 
correspond to low energy excitations in the vicinity of the minima (and therefore 
will be unimportant for local properties like Green's functions), but will be essential 
for smoothing the singularity at the cusps. In the example just discussed one would 
have to include in the game at least A^ fields. That is how A^ enters the effective 
potential as a hidden parameter (besides the overall A^^ factor in Eq. (^)). Then 
the typical value of relevant fields Mj inside the wall solution can be 0(A^~^), each 
field contributes to F at the level 0{N~^), but there are A^ relevant fields, and the 
value of F ~ 1/A^. Note that the wall width is L ~ N^. Then the volume energy 
density inside the wall E ^ N and, correspondingly, the wall tension e ^ N. 

The fact that the volume energy E is 0{N) inside the BPS wall connecting two 
neighboring vacua, say with Tr(AA) = exp{27iik/N) where k = and 1, is seen 
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in the microscopic theory (|^) per se. Indeed, for the BPS wall 




Since the volume energy density E ^ N G^j^G^^,, and AA^ in the neighboring chirally 
asymmetric vacua is 0(1), we conclude that E scales as A^. 

How do we learn about the A^-sector structure of the effective Lagrangian emerg- 
ing in SUSY gluodynamics? If the gauge group is SU{N), the theory has a discrete 
Z2N chiral symmetry, which is spontaneously broken down to Z2 in some of the 
vacua. This means that the effective Lagrangian for mesons/glueballs must have 
at least degenerate minima which differ from each other only in the value of the 
phase (J) of the order parameter TrA^ (the latter is an interpolating field for one of 
the lighter mesons). The minima lie at = 0, 27r/A^, 47r/A^ and so on. Then, clearly, 
we must have a much more rapid variation of F as a function of 0, than one would 
naively infer. Naively, since there is no explicit N dependence in E, one would say 
that if the first minimum lies at = 0, the second one should be at ~ 1. The 
only way out is either to have an A^-sector structure (within the construction that 
includes only fixed, A^ independent number of fields in the effective Lagrangian), or 
to build a Lagrangian on a minimal set of A^ fields. In both cases a hidden parameter 
A^ appears. It does not affect the value of the derivatives of E at the minima, which 
are all of order 1. 

In Sect. 3 we will consider the amended Veneziano-Yankielowicz effective La- 
grangian and will see that it indeed has the required structure. In a somewhat 
simplified picture, we can understand how it appears by considering the dependence 
of the vacuum energy in the Yang-Mills theory on the vacuum angle ^9. As is well 
known form the consideration of the Ward identities |T^, this dependence has a 
"wrong" periodicity. That is, naively the energy is periodic in i} with the period 
2ttN rather than 27r, 



The correct periodicity of the physical quantities is restored in the following way. 
The ground state at < vr, becomes an excited state at > ir. At = ir there are 
two degenerate states. At this point due to the level crossing the vacuum energy 
has a cusp so that 



Now, the interaction of the phase with the gluonic degrees of freedom is the same 
as of the rescaled angle i)/N. The effective potential for is, therefore, roughly 




(7) 




(8) 



i?vac OC A^^ —1 + COS — 



2TTk 



N 



) 



{2k - l)n 
N 



< < 



{2k + l)n 
N 



(9) 
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This potential indeed has the form of Eq. (^. Moreover, the derivatives of the 
function F at all the minima are 0(1). Nevertheless, it has N minima at = 2TTk/N . 
The value of F on interpolating trajectories varies from zero at the minima to 
0{N~'^) at the cusp. A naive estimate of the wall tension would therefore give a 
value much smaller Qthan 0{N^). 

Our next remark concerns a field-theoretical understanding of a confining string 
which ends on a domain wall. A simple example of such a situation is the wall 
that separates the confining phase in a gauge theory from a nonconfining one. The 
type I wall in SUSY gluodynamics is precisely of this kind. Recall that the chirally 
symmetric vacuum at (A^) = sustains massless excitations. It was argued in 
that this phase is in fact conformally invariant and can be thought of as a kind 
of a non-Abelian Coulomb phase. Now consider a wall that separates a confining 
phase from a Coulomb phase. A probe charge placed in the confining phase is a 
source of the electric flux which travels in a flux tube - the confining string. On the 
other side of the wall, however, it is energetically favorable for the flux to spread 
out into a Coulomb field. So an observer in the Coulomb phase will not see a string 
but, rather, a point charge (in fact, twice as big in magnitude as the original probe 
charge, since the electric flux will spread in half the space) sitting on the domain 
wall. One is not used to thinking about a phase boundary between the Coulomb 
and confining phases, since usually the two are not degenerate in energy. It is the 
peculiar feature of supersymmetric theories that two physically completely distinct 
phases are degenerate. 

It is very easy to find an easy-to-handle example of a domain wall separating the 
confining and the Abelian Coulomb phases by introducing some extra fields in SUSY 
gluodynamics. Start from the Lagrangian (§), and add one chiral matter superfield 
$ in the adjoint representation of the gauge group, with the superpotential 

>V = mTr$2 + Af"^(Tr$2)2 (10) 

(the gauge group SU{2) is assumed). The second term in the superpotential is non- 
renormalizable. One can think of it as a result of integrating out some heavy degrees 
of freedom, so that at a large scale M we return back to a renormalizable theory. It 
is assumed that m <^ A, but mM S> A^. 

If = the theory we deal with is nothing but a softly broken N = 2 



model studied by Seiberg and Witten |16|. In the Seiberg-Witten vacua Tr$ ~ A 



monopole condensation takes place, and due to the dual Meissner effect the probe 
electric charges placed in one of these vacua will form flux tubes. The presence of a 



very weak additional interaction M~^(Tr$^)^ not considered in Ref. [16| does not 
affect the picture obtained there, since this term can be viewed as an arbitrary small 
perturbation if the theory resides in one of the Seiberg-Witten vacua. 

However, at large values of Tr$^ the term M^^(Tr$^)^ leads to a drastic re- 
structuring - no matter how small is there appears a new vacuum state at 

^In fact, a naive (and wrong) estimate would give 0{N'^). Actually, a much larger contribution, 
0{N^), resides in the cusp, see Sect. 3. 
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Tr<l>2 ~ mM. In this vacuum the gauge symmetry is broken down to [/(I) by a very 
large vacuum expectation value of the $ field, the monopoles are very heavy, and 
the theory is obviously in the weakly coupled Coulomb phase. Supersymmetry guar- 
antees that the vacuum energy densities in both phases vanish: the two phases are 
degenerate. Under the circumstances a domain wall separating the weakly coupled 
Coulomb phase and the strongly coupled confining phase (one of the Seiberg-Witten 
vacua) must exist, with the wall tension ~ m?M. If the confining phase is to the 
left of the wall, and we put there a probe electric charge, a fiux tube going towards 
the wall develops; the chromoelectric fiux is clearly diffused to the right of the wall. 

Another example of a wall that serves as a sink of the chromoelectric flux is a 
situation when the Coulomb phase exists not in half the space (say, to the left of 
the wall) but only inside the wall. For example if one considers a wall in SUSY 
gluodynamics that separates the phases TrA^ = and TrA^ = — A^, it is very 
likely that the order parameter A^ will vanish inside the wall. Both phases then 
are conflning but the wall itself is "made" of the Coulomb vacuum. In such a case 
the flux tube that originates in one of the phases will not penetrate into another 
but, most likely, the flux will spread out in transverse directions inside the wall in a 
two-dimensional Coulomb fleld. Energetically this is preferable, since the energy of 
the two-dimensional Coulomb fleld depends only logarithmically on the size of the 
system, while the energy of the string that penetrates into the other phase is linear. 

Note that in this latter scenario the degeneracy between the Coulomb and con- 
flning phases is not necessary. The picture can be dynamically realized both in the 
super symmetric and nonsupersymmetric contexts. An illustrative nonsupersymmet- 
ric model, where inside the wall the theory is in the Abelian Coulomb phase while 
outside it is in the conflning phase, was presented in Ref. It would be instructive 
to exploit this model for a more quantitative analysis of a flux tube coming from the 
conflning phase and diffusing itself inside the wall (i.e. in the Coulomb phase). A 
(semi) quantitative analysis seems possible since at least inside the wall the theory 
01 is in the weak coupling regime. 

Finally, if in the previous examples we consider the Higgs phase instead of the 
Coulomb phase (either to the left of the wall, or inside it), the chromoelectric flux will 
still disappear in the wall. In this case it is even more trivial, since the chromoelectric 
flux is not conserved in the Higgs phase, and it will be screened by the Higgs phase 
vacuum either on the left side or inside the wall. We would like to argue that 
the type II wall, considered in Ref. 0, is, in fact, an example of this kind, in a 
certain sense. Indeed, consider the type II wall. Let us say, on the left there lies 
a phase with the condensate of monopoles, while on the right with the condensate 
of dyons 0. Let us imagine a probe electric charge to the left of the wall. Since 
the dyons are electrically charged, their condensate acts like a Higgs vacuum, in 
the sense that it can be easily polarized to completely screen the electric flux that 
might enter the dyon condensate through the wall. Of course, since the dyons are 

■^Both, the monopoles and dyons in SUSY gluodynamics are to be understood in the same sense 
as those in the 't Hooft construction . 
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also magnetically charged, any such polarization of their condensate will lead to 
appearance of net magnetic charge to the right of the wall. However, the magnetic 
flux tube emanating from this induced magnetic charge can be directed towards 
the domain wall. In that case it will be screened on the other side of the wall by 
polarization of the condensate of the monopoles. In other words it is plausible that 
the dyonic condensate to the right of the wall will be polarized to screen the electric 
charge while the monopole condensate to the left of the wall will be polarized to 
compensate for the excess of the induced magnetic charge. As a result the confining 
electric string will terminate on the wall. 

Note, that this picture is somewhat different from the one advocated in Ref . , 
where it was suggested that a bound state of a monopole and a dyon appears on the 
wall. It is difficult to talk about monopoles and dyons forming a bound state, since 
they do not exist as free particles on either side of the wall, because both vacua have 
nonvanishing condensates. 



2.2 Toy model — axion wall 

The type II wall can be thought of as carrying the quark quantum numbers in the 
presence of a QCD string, in the sense that it can screen a fundamental charge or 
other charges which are nontrivially transformed under the center Z^v, in the theory 
where all dynamical fields are invariant with respect to Z^^. Surprising as it is, one 
can trace the very same phenomenon in simple Abelian models. Although the paral- 
lel is not perfect, a simple Abelian example may be useful for deeper understanding 
of this general phenomenon. 

The problem we keep in mind is the axion wall in the presence of a monopole 
||18|| . For simplicity we will consider the SU{2) case (the Georgi-Glashow model). 
The SU{2) symmetry is spontaneously broken by the vacuum expectation value of a 
Higgs field down to U{1) giving rise to the 't Hooft-Polyakov monopoles |jl9|]. After 



the breaking, the fields in the adjoint representation have the U{1) charges ±1, while 
those in the fundamental representation have charges ±1/2. 

Let us recall some facts about the monopoles in the presence of the term. The 
Lagrangian of the Georgi-Glashow model is 

C = -^FrF;:^ + ^,FrF^^ + Lh($) , (11) 

where e is the gauge coupling constant and the last term is the scalar Lagrangian 
for the Higgs field in the adjoint representation of the SU{2) group. 

It was shown by Witten that if ^ 0, a monopole becomes a dyon with the 
electric charge 



87r2 



/i + ne , (12) 
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where /i is the magnetic charge, 



4vr , , 

/x=— . 13 

e 

When d changes from to 27r one gets n ^ n + 1. 

In the theory (|TT]) "Q is constant, given once and forever. However, if the axions 
are added in the theory, then, effectively, "d is substituted by the axion field which 
can vary in space-time. The axion field can be introduced through a phantom-axion 



construction i.e. we add an SU (2)-singlet Higgs field coupled to a doublet quark 
field. In the limit when the expectation value of the 5'[/(2)-singlet Higgs field tends 
to infinity, the quark becomes infinitely massive and disappears from the spectrum, 
and so does the modulus of the singlet Higgs field. Its phase becomes an axion field 
a{x). 

The axion Lagrangian is 

l^a = ^-fd^ad^a -K^[l- cos(a - d)] , (14) 

where the parameter K"^ is connected with the vacuum susceptibility, it is exponen- 
tially small in the model at hand, ~ exp(— Svr^/e^). Moreover, fa is (a very large) 
expectation value of the singlet Higgs. In this limit the only other axion interaction 
to be taken into account is its coupling to FF. The d term in Eq. (|lT]) becomes 

^^ = ^^i^^"^-- ^^^^ 

One can now set ■(9 = (and we will do this hereafter). 

The vacua at a = and a = 2'k are physically identical. Correspondingly, axion 
domain walls exist in this system ||22| interpolating between a = and a = 2tt. 
Assume a wall lies in the xy plane, so that the axion profile depends only on z. The 
wall solution centered at origin is 

a{z) = 2-71 — 2 arccostanh(m2;) (16) 

where m is the axion mass, m = K/fa, and the width of the wall is of order m~^. 
Of course, it is assumed that m <^ M where M is the monopole mass. 

Start from the monopole with electric charge zero to the left of the wall, and 
let it adiabatically propagate through the wall. Effectively, adiabatically changes 
from to 27r. To the right of the wall the monopole becomes a dyon with electric 
charge 1. Inside the wall, the electric charge of the monopole gradually increases. 

Thus, one gets an apparent nonconservation of the electric charge of the monopole. 
However, since U{1) is unbroken , the total electric charge must be conserved. The 
question is where is the missing electric charge. 

As was shown by Sikivie the monopoles actually induce electric charges 



±1/2 on the wall. When the monopole is far to the left of the wall, it is neutral, but 
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the charge induced on the wall is +1/2. When the monopole is far to the right of the 
wall, the induced charge will be —1/2, so the total charge is conserved, 1/2 = 1 — 1/2. 

To see that this is indeed the case we observe that the extra term (|r^) in the 
Lagrangian is immediately translated into an additional piece in the definition of 
the electromagnetic current 



■axion _ "-^ 
J u 



5A^ Stt 



^^^Kl) (17) 



where it is assumed that the expectation value of the triplet Higgs field is aligned 
along the third direction, so that F^^ is the photon field strength tensor. Note that 
jaxion -g automatically conserved, Q^j^^"^ = 0. The corresponding contribution in 
the electric charge consists of two parts, 

Q = J d'xjr'' = ^ {Fld,<^ + ad^Fl) . (18) 

Let us assume that the distance between the wall center R and the monopole is much 
larger than m~^. For such distant monopoles the physical meaning of each term in 
Eq. (|18D is transparent. The second term vanishes everywhere except the point 
where the monopole sits. Thus, it gives the electric charge of the monopole/ dyon. 
If the latter sits to the left of the wall, where a = 0, the monopole electric charge 
vanishes. To the right of the wall a = 27r, and 

^Qmonopole J d 3^V-Bnionopole 1 • 



The first term in Eq. (^8|) is obviously saturated inside the wall; it describes the 
electric charge induced on the axion wall in the presence of a distant monopole. 
The induced charge is equal to the flux of the monopole magnetic field through the 
plane of the wall times Aa/(87r^). Since this flux is 1/2 of the flux through the large 
sphere (= e/i/2 = 27r), the induced charge on the wall is obviously equal to ±1/2 
depending on whether the monopole is on the left or on the right of the wall. 

Thus, the picture is in complete agreement with the conservation of the total 
electric charge. 

This picture can be readily generalized for SU{N). Then there are A^ — 1 different 
monopoles corresponding to — 1 (= rank for SU{N)) Abelian U{1) factors in 
the Cartan subalgebra of SU{N). Repeating the same analysis, one can see that 
fundamental and antifundamental representations [A^] and [N] are induced on the 
domain wall (n, n+1) and (n + l,n) respectively. Here n = 0, 1, ...A^ — 1 corresponds 
to a = 27m. One can see that taking two domain walls (n, n + 1) and (n + 1, n + 2) 
one can get all representations corresponding to the product of two fundamental 
representations [A^] x [A^]. For m consecutive walls we have [A^] x [A^]... x [A^]. If one 
takes m = A^ — 1 than one can get antifundamental representation, in full agreement 
with the fact that the A^ - 1 walls (0, 1), (1, 2), ...(A^ - 2, A^ - 1) make together a 
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wall (0, — 1), which is equivalent to the wall {N, N — 1). But the last one, as any 
(n + l,n) wall, corresponds to antifundamental representation [N]. 

Returning to our original problem, SUSY gluodynamics, we note that the phase 
of in a sense plays a role analogous to the axion field. It could be interesting to 
pursue the analogy between the Abelian toy model and SUSY gluodynamics further. 

In summary, in supersymmetric theories which have degenerate vacua with very 
different physical properties, the fact that the confining string can end on a domain 
wall is quite natural. Actually, the wall does not have to be BPS saturated to serve 
as a sink for the chromoelectric flux carried by the string. 

Regardless, it is still an interesting question whether all possible BPS saturated 
walls are dynamically realized in SUSY gluodynamics. In the next section we will 
attempt to address this question. Although we will not be able to give a posi- 
tive proof, we will show that the straightforward search in the framework of the 
Veneziano-Yankielowicz Lagrangians for the walls connecting neighboring chirally 
asymmetric vacua is in general a dangerous endeavor. As we shall see the cusp 
structure of these Lagrangians makes it impossible to decide this question without 
additional nontrivial dynamical information. We will also present a toy version of 
the underlying phenomenon. 



3 Glued Potentials 

3.1 counting and paradoxes of the wall-building in the 
Veneziano-Yankielowicz Lagrangians 

First, we briefly remind the relevant formalism. The effective Lagrangian for SUSY 
gluodynamics was written down a long time ago and then amended recently |^ 
to properly incorporate the non-anomalous symmetry. 

We will write down the Lagrangian realizing the anomalous Ward identities in 
terms of the chiral superfleld 

S= TtW\ (19) 

namely. 



^ j d^ONS i\n ^ + 2Tin^ + d^ONS (\n ^ - 27rm'j , (20) 



3 7 V ^ / 

where a is a numerical parameter, 

a = eA^^e^''/^ 
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A is the scale parameter, a positive number of dimension of mass which we will set 
equal to unity in the following. Please, note the factors in Eqs. (p!9D and (pO|). 

An important element in the Lagrangian ( PPD is an integer-valued Lagrange mul- 
tiplier n. In calculating the partition function and all correlation functions the sum 
over n is implied. The variable n takes only integer values and is not a local field. 
It does not depend on the space-time coordinates and, therefore, integration over it 
imposes a global constraint on the topological charge. It is easy to see that (after 
the Euclidean rotation) the constraint takes the form 

^ = j d'xG%G% = Z . (21) 



While the F term in Eq. ( pOD is unambiguously fixed, the D term is not specified 
by the anomalous Ward identities. We have chosen it in the simplest possible form, 
with the numerical coefficient which gives the correct large A^ counting. 

The extra term added to the Lagrangian is clearly supersymmetric and is also 
invariant under all global symmetries of the original theory. The single-valuedness 
of the scalar potential and the Z^r invariance which were missing in the original 
Veneziano- Yankielowicz effective Lagrangian are restored 0. The chiral phase rota- 
tion by the angle 2TTk/N with integer k just leads to the shift of n by units. Since 
n is summed over in the functional integral, the resulting Lagrangian for S is indeed 
Ztv invariant. 

The constraint on [5— following from the Lagrangian ( pO|) results in a peculiar 
form of the scalar potential. The expression for the scalar potential is given in Ref. 
i- 

Eliminating, as usual the F component of S with the help of classical equations 
of motion at fixed n, the effective potential can be written as 



f/(0) = -r-Mn 



^exp{-4Ar2v(0*0)2/3[in2 |0| + (« + vm)^]} 



f22) 



Here V is the total space-time volume of the system, cf) is the lowest component of 
the superfield 5*, and a = Arg(0). In the limit V ^ oo only one term in the sum 
over n contributes for every value of a; which particular term depends on the value 
of a. Thus, for —n/N < a < tt/N the only contribution comes from n = 0. In this 
sector the scalar potential is 

U{(f)) = 4A^2(0*0)2/3 In In 0* , (23) 

In general we have 

f/(0) = 4Ar2(0*0)2/3ln(0e-'2WA^)ln(0*e*2WA^) (24) 

^The explicit invariance here is Zpf rather than the complete Z2N of the original SUSY gluody- 
namics, since we have chosen to write our effective Lagrangian for the superfield which is invariant 
under A — s- —A. 
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(2n - 1)71 , (2n + l)n 

at ^ <arg</)< ^ . 

In other words, the complex plane is divided into sectors. The scalar po- 
tential in the n-th sector is just that in the first sector rotated by —27t/N. The 
scalar potential itself is continuous, but its first derivative in the angular direction 
experiences a jump at aigcf) = {2n + l)n/N. The scalar potential is "glued" out of iV 
pieces. The Zjy symmetry is explicit in this expression. It is quite obvious that the 
problem at hand has + 1 supersymmetric minima - N minima at = e*^'^"/^ , cor- 
responding to a non- vanishing value of the gluino condensate (spontaneously broken 
discrete chiral symmetry), and a minimum at = (unbroken chiral symmetry). 

Including the kinetic term of the field S, as it appears in Eq. (|20|), leads to the 
following effective Lagrangian: 

L = N^d,<p'/%<P*'/' + Um. (25) 

We now ask ourselves whether this Lagrangian can be used to find an explicit 
wall solutions. In fact, the solution for the type I wall has been considered in detail 
in [0 and was found to exist and to be BPS-saturated. The situation with the type 
II walls is more complicated. Note that any field configuration that interpolates 
between the two vacua at = 1 and = e*^'^/^ has to go through a point where the 
phase of the field is vr/A^. At this point the scalar potential has a cusp, and one 
has to be very careful in treating such configurations. 

As an illustration, let us forget for a while about possible complications and 
estimate the tension of the type II wall using the Lagrangian (pSj). The potential 
and kinetic energies should contribute to the tension of the wall roughly the same 
amount, so we concentrate on the kinetic term. The variation of the field inside 
the wall at large A^ is A0 ~ 0{1/N). The mass m of the field is independent 
of A^ with our choice of the coefficient of the kinetic termQ. The width of the wall 
obviously is of order L ~ 1/m ~ 0(1). The kinetic energy contribution is, therefore, 

Ekin ~ N\A<f)f/L' ~ 0(1) , e ~ Ekin^ ~ 0(1) . (26) 

The same estimate is obtained if we consider the potential energy contribution. 

Surprisingly, this is by far lower than the BPS bound on the wall energy in the 
original theory, Eq. (^, which gives e ~ 0{N). At first sight this might seem to be 
an arithmetic paradox even in the framework of the effective theory per se, viewed 
as a generalized Wess-Zumino model. In the generalized Wess-Zumino models with 



a superpotential yV{S) the BPS bound is (0, §], see also p4| ) 



£>2|W(5i)-W(S'o)|, (27) 

where 6*0,1 are the values of the field in two vacua between which the given wall 
interpolates. Taken at its face value that would give a bound of e ~ 0{N) for the 
type II wall we are considering. 

^This is, in fact, how the meson masses should behave at large N and this is the reason of 
choosing the coefficient A^-^ in front of the kinetic term in Eq. (|20|). 
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In fact, there is no arithmetic paradox here. The BPS bound (|27| ) on the wall 
energy in the generalized Wess-Zumino models assumes that the superpotential is 
smooth. In the effective theory (PD|), for the walls that cross the cusp (type II), it 
experiences a jump, 

A>Va=W(5a+)->V(5'a-)^0. 

(Here S'a is the point where the wall crosses the discontinuity line, and >V(S'a-) and 
>V(S'a+) are the values of the superpotential below and above this line, respectively.) 
The superpotential which is obtained from Eq. (|20|), after summation over tt,, has 



a phase discontinuity along the same lines where the effective scalar potential 
has a cusp. Accounting for the jump, modifies the bound, and, instead of Eq. 
we, therefore, have 

2 \W{Si) - W{So) - AWaI . (28) 

Due to the discontinuity AWa, the expression in Eq. (|28|) is 0(1), in full accord 
with Eq. (^. 

Nevertheless, even though there is no paradox at the level of the effective theory, 
clearly the effective potential grossly misrepresents the tension of the type II wall. 
This is surprising since the field (p changes slowly inside the wall, and normally 
one would think that the effective potential should properly describe slowly varying 
configurations. It is clear that this failure is intimately connected with the cusp 
structure of the effective potential. Our aim now is to understand what is the 
physical origin of the cusp structure. To warm up we consider first a very simple 
(non-supersymmetric) model, leading to a similar structure, and then move on to a 
more general picture of the phenomenon in supersymmetry. 



3.2 Glued effective potential in a simple model 

To understand how an effective potential with cusps can appear from a smooth 
potential of the original theory it is best to consider an explicit example. Let us 
take a (non-supersymmetric) theory of two scalar fields with the potential (Fig. 1) 

U{<P, X) = ^ i<P' - vr + ^ (X' - - '-'<PX (29) 

where and x ci-re real fields and A and ( are the coupling constants. The coupling 
constant k is taken to be real. The mass of the (p quantum = 4A?7^, while that 
of the X quantum = 4:(fi'^. Let us assume that the field is much heavier than 
the field x, M ^ m. For technical simplicity we will also assume that k <^ m, M, 
and the expectation values 1], fi are of the same order of magnitude, although this is 
not crucial. 

The theory has two symmetry breaking minima. When both coupling constants 
A, ( are large, those are 

{<P = V, X = /^} and {0 = -7], x = -/^l ■ (30) 
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Figure 1: The scalar potential in the model considered in Sect. 3.2. 



Let us now derive an effective potential for the field x by integrating out (p. To 
calculate the effective potential in the leading adiabatic approximation we fix the 
value of X ^^id solve for cf) in this background. Note that for a fixed and not too 
large value of x the potential for </> has two local minima. Generically, the two 
local minima are nondegenerate. For x < the state = —r] has lower energy 0, 
while for x > the global minimum is at = r^. At x = both local minima 
become degenerate. At this point there is a discontinuous change of the vacuum in 
the heavy sector. As a result the effective potential develops a cusp of precisely the 
same nature as discussed in the previous section, 

U,^ = -K^r^\x\ + \{X^-I^f. (31) 

Thus, although the underlying theory is perfectly smooth, the effective theory has 
two sectors and a cusp at x = (Fig. 2). 

We can again pose the question about the existence of the wall, and the cal- 
culability of its energy density from the effective Lagrangian. Clearly there exists 
a solution of the original equations of motion, stemming from Eq. (p9D , with the 
wall boundary conditions, i.e. interpolating between the two vacua of Eq. (^). 
It is equally clear that the estimate of its energy from the effective Lagrangian 
will be generically incorrect. The obvious reason is that the effective Lagrangian 
is completely independent of the larger mass scale M . This is natural since it was 
calculated in the leading adiabatic approximation, i.e. in the limit M oo. On the 
other hand, to produce the wall one has to excite the heavy field 0, which jumps 
from —7] to +7] inside the wall profile. This costs energy proportional to M, so the 
wall energy density in this theory must be proportional to M. The wall tension in 
the present example can be calculated directly from the "fundamental" Lagrangian 

^We neglect here small corrections of order /i/77 and K/rj to the values of at the minima. 
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(p9|) without appealing to the effective Lagrangian (^TJ). Roughly it behaves as 

e = xMrf + ym/i^ (32) 

where x and y are numbers of order one. If the expectation values of the heavy and 
light fields are of the same order, ~ /i, the bulk of the wall energy is contributed 
by the heavy modes. In that case the wall energy can not be obtained from the 
effective potential for x- O^is can say that a dominant part of the wall tension is 
associated with the cusp. 

Physically the picture of what is happening is very simple. The field x is light 
and therefore changes slowly inside the wall, on the scale 1/m. The heavy field 
follows this change adiabatically almost everywhere in space, except for the region 
where x = 0. In this region, within a distance of order 1/M, the value of changes 
from —7] to 7]. The big contribution to the wall energy density, proportional to M 
comes precisely from this small region in space in which the field x sits on the cusp 
of the effective potential. 

If we use the effective potential to calculate the wall energy, the result will be of 
order m/i^, since this is indeed the contribution of the light field. The contribution 
of the heavy fields can be thought of as an extra contribution of the cusp in the light 
field effective potential. 

There are several lessons we want to draw from this toy model. First, the physical 
reason for the appearance of the cusp in the effective potential (glued potential) at 
some value Xo of the field is that at this particular value of the light field, the system 
of heavy fields has two (or more) degenerate ground states. In general, when the 
value of the light field changes continuously, the heavy fields follow adiabatically 
after this change. However when the light field passes through xoi there is a level 
crossing in the heavy system and the properties of the heavy field vacuum change 
discontinuously. This "first order phase transition" leads to discontinuity in the first 
derivative of the effective potential, and, therefore, a cusp. 

The second lesson, is that the wall configuration that connects the points on 
different sides of the cusp necessarily involves excitation of heavy modes. This is so 
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since close to the point Xo inside the wall the heavy modes must rearrange in order 
to make the transition between the two degenerate vacuum states. 

Finally, the description of the wall with the help of naive glued potential is not 
valid if the bulk of the wall energy comes from the change of the heavy fields. The 
problem stems from the fact that the effective potential is calculated in expansion 
in powers of 1/M, which is the adiabatic approximation. As we already stressed, 
the adiabatic approximation breaks down inside the wall due to the level crossing. 
On one hand, this is a signal of possible appearance of terms proportional M in 
the effective action. On the other hand, this is precisely the situation in which 
a nontrivial topological Berry phase should appear. A more careful calculation of 
effective action should reveal the presence of the terms of the type 

M[(/)(oo) - (/)(-oo)]^ (33) 

This topological term does not affect the vacuum sector, but adds the missing large 
piece, cusp contribution to the energy of the wall. 

Returning to SUSY gluodynamics it is now clear why the original BPS bound 
is so badly violated by a wall configurations in the effective theory. The reason is 
that the effective theory misses a large cusp contribution to the energy which comes 
from the heavy modes not appearing in the effective Lagrangian that are excited 
inside the wall. We conclude, therefore, that the Veneziano-Yankielowicz effective 
Lagrangian, as it stands, can not be used to calculate the wall energy; without 
additional information we can not say whether or not the BPS saturated type II 
walls exist in SUSY gluodynamics. 

In the next section we would like to give a detailed example of how this situa- 
tion arises in a supersymmetric theory. The model we will consider has the same 
symmetries as SUSY gluodynamics and an effective potential of the Veneziano- 
Yankielowicz type. We will see in detail how the cusp structure of the effective 
Lagrangian appears when integrating out the heavy superfields and will be able to 
trace exactly the missing piece of the wall energy. 



3.3 Supersymmetric model with sector superpotential 

The model of the previous subsection was only intended for explaining how cusps 
arise in the effective potential. We now want to consider a model which captures 
more features inherent to SUSY gluodynamics. Consider a generalized Wess-Zumino 
model of two scalar chiral superfields $ and X, with the superpotential 

W = - b<^X + ^X' , (34) 

where the coefficients a, b and c are real positive numbers. This model obviously has 
a Zjv symmetry, under which 

^^^^i2nn/N^ X^Xe*^""/^, ^^^e*^'^"/^. (35) 
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The global minima of the energy are determined from the equations 



^ = a*"«-W = 0, ^ = -(,* + cA-=0. (36) 

09 oX 

These equations have one symmetric solution 

$ = X = (37) 
and N solutions which spontaneously break the symmetry 

$n= - e^^-"/^ X„ = -<|.„. (38) 
\acj c 

Choosing 6 ^ c, we find that the field $ near the asymmetric vacua is very heavy; 
its mass 

M=(X+1)-, (39) 
c 

while the field X is light, with the mass m = c 0. 

The effective potential for the field X is obtained by eliminating the heavy field 
$ by virtue 9VV/(9$ = 0, see the first equation in (|36D. The condition (9W/9$ = 
has N + 1 solutions 

(h r 27rm 1 , , 

MX)=[-X) exp{ — }. (40) 

The solution has to be substituted in the superpotential, 

W[$,X] ^ W[$*(X),X] = 




l/(iV+l) 



^(N+2)/(,N+l) 



( 27rin \ c 



(41) 



This is not the end of the story, however. The effective superpotential is obtained by 
choosing for every value of X the solution that gives a minimal energy. The energy 
as a function of the lowest component x of fhe superfield X for each branch is 




exp 



2nin 
N+1 



+ cx 



(42) 



^Strictly speaking, the mass matrix is non-diagonal. There is a smah admixture of X in the 
heavy diagonal combination, and a small admixture of $ in the light diagonal combination. These 
admixtures are 0{cb~^N~^), and the absolute shift in the mass eigenvalues is 0{cN~^). This shift, 
as well as the mixing mentioned above, can be neglected in the limit iV 3> 1 and b ^ c. Moreover, 
if 3> 1, it is not necessary to require b ^ c. Even at 6 ~ c the hierarchy of masses takes place, 
M ^ m, and the off-diagonal elements of the mass matrix are negligible. 
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Clearly, the energy is minimal for the branch for which Argx is closest to 



27cn 1 , 

-Arg X ■ 



N+1 N+l 

For instance, if x is real and positive, = 0, if Argx = 27t/N then n = 1, and so 
on. At Argx = tt/N both branches, with n = and n = 1, have the same energies. 
The resulting effective superpotential has sectors and is "glued" along A^ rays. 

We. = lx'-^^b , (43) 

TT TT 

at < ArgY < — , 

at — < ArgY < — , 

and so on (see Eq. (|41|)). The discontinuities in the effective superpotential (or, 
equivalently, the cusps in the effective potential) occur along the rays 

(2A; + l)7r 
Argx = j;^ 

where two branches in Eq. (|40|) with n = k and n = k + 1 are degenerate in energy, 
see Eq. (^2]). For instance, at Argx = n/N 



>Veff(XA±) = e 



2ni/N 



C|XP N+l ( h\ |(7V+2)/(iV+l) " 



Now we are ready to address the issue of the domain walls. Consider a domain 
wall that connects two adjacent asymmetric vacua. The BPS bound on its tension 
is 

2A^ 62 /^2\ # _ .27r\ 



^>2|>V(».,.^o)-m*..A-0l = ^-^-j sm^-J. (45) 

At large A^ this is of order 1/A^. On the other hand, a naive estimate based on the 
effective potential for the light field would yield 

^2 / ^2 \ :! ]^ 

,.„|A,f _ -^M-\ (46) 
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Just like in SUSY gluodynamics the two expressions are incompatible. We know 
already that the reason is that the adiabatic approximation used to derive the effec- 
tive Lagrangian breaks down at the cusp. For configurations crossing the cusp an 
extra "topological" term has to be added to the effective Lagrangian, as discussed in 
the previous section. Let us first estimate a part of the tension associated with the 
cusp, and the corresponding restructuring of the heavy field (f). A straightforward 
estimate analogous to that of Eq. (^61) is 




£^ ~ M|A0|^ ~ - - T-r^ (47) 



i.e. we reproduce the order of magnitude of the BPS bound (|45|). 

As a matter of fact, in the present case one can do better than that. To this end 
we note, that the only point where the adiabatic approximation breaks down is at 
the cusp. In other words, almost everywhere throughout the space the heavy field (f) 
does indeed follow the change of x adiabatically. Only at the point inside the wall, 
where 

X = pe-/^ , 

the value of changes rapidly. This change, of course, does not happen abruptly, 
but rather on the scale of the inverse mass of the field 0. The field x remains 
constant throughout the region of space where the rapid variation of (p takes place. 
The profile of in this region as well as the energy associated with this variation 
can be calculated by considering the original theory at a frozen cusp value of x- 

lfx = P exp {iir / N) , the wall profile of the field is determined from the following 
superpotential 

At X = pexp {m/N) the two branches of Eq. pO|) are degenerate. The wall under 
consideration interpolates between 




5*0,1 = \-p\ exp 



, 1 1 

iix \ — ±. 



N N +1 



(49) 



where the upper and lower signs correspond to the final and initial points, respec- 
tively. These points are two degenerate minima of the potential stemming from Eq. 
(^). The wall in question is BPS saturated. This is because the BPS equation in 
this case is a pair of the first order equations for two real fields (real and imaginary 
parts of 0), 

dW f in 2'Ki\ , , 
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which possess one conserved quantity (see Ref. Q for details) 



Im 



/ivr 


27riY 




1^). 



0. (51) 



Using Eq. ( ^iD one can always eliminate one real field, getting in this way a trajec- 
tory in the plane {Re0 , Im 0} that connects 0*o,i- The trajectory depends on one 
real variable. The resulting one first-order equation for one real function always has 
a solution. 

We conclude that the tension of the wall (the cusp term) is given by the BPS 
bound for the theory with the superpotential 



This is precisely the term that has to be added to the effective potential of x 
in order to be able to calculate the wall tension properly. Equation (|5^ can be 
generalized to a wall configuration which connects any two vacua, and not neces- 
sarily the two adjacent ones. The extra term in this case would be just the sum of 
contributions of all cusps crossed by the wall. 

Note that although the component is built above through the BPS saturation, 
the full {(f) , x} wall in the original theory is not BPS-saturated, at least, in some 
range of parameters. (We mean the type II wall connecting, say, two neighboring 
asymmetric vacua with n = and n = 1 in Eq. (pSf).) 

To see that this is indeed the case consider the model (^) with a = b = c = 1 
and ^ 1. (For these values of parameters the field $ is still much heavier than 
X, see the footnote following Eq. (p9D). Assume that the wall is BPS saturated. 
Then 

d.^ = (0^+^ - x) exp (f + ^) , d.x = (-0 + X) exp (f + ^) • (53) 
As a consequence. 



5.0 + d^x = [r^' - 0j exp + — j , (54) 

where the right-hand side contains no dependence on x- We can take advantage of 
this fact. Multiply both sides by ^20 and integrate over z from — oo to -|-oo. Then 



I-oo^^"^ + d,x)dz<pdz = exp (^y + 



IN +2 _ 1 12 



n oo 



N + 2 



(55) 



where ^bps is the tension defined on the right-hand side of Eq. (^q ) and the values 
of the field in the vacua, Eq. (p8l), are substituted. Equation (^) implies 

{dzX){dJ)dz = -'-^-C, (56) 
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where Ci is a real positive number. On the other hand, 

-oo J —oo 



{d,x){d,^)dz + }i.c. 



(57) 



where C2 is a real positive number. The first term is obviously equal to eBPs/2, 
while the second term can be read off from Eq. (|56|). In this way we obtain 



Co 



^BPS 



+ 2 



(5J 



or C2 + 2Ci = — £:bps/2- This relation is obviously inconsistent, which proves that 
we cannot built a BPS saturated wall in the problem at hand. 

Thus, the lesson to be drawn is as follows. Consideration of the effective low- 
energy theory, by itself, yields information on the vacuum structure. The fact of 
existence of the walls can be unambiguously inferred from this information. But 
neither the nature of the wall (BPS versus non-BPS), nor its tension can be prop- 
erly found from the analysis of the effective low-energy theory if the corresponding 
potential is glued from distinct sectors, and the wall in question crosses the cusps. 



3.4 Elements of the general theory 

Given an effective low-energy theory, obtained after integrating out all heavy fields, 
with a discrete set of degenerate vacua (as it is typical for supersymmetry) , the 
question we ask is: can one infer from this low-energy theory the existence of the 
BPS walls interpolating between the distinct vacua? Under what circumstances the 
BPS wall seen in the effective theory is a reflection of the wall in the full theory? 
And vice versa, if we see no BPS walls in the effective low-energy theory does it 
mean there are no such walls in the full theory? 

The full general theory is not yet worked out, and the answers to these questions 
in the generic situation are not known so far. In this section we will present some 
illustrative considerations which are valid in the simplest possible setting: the gen- 
eralized Wess-Zumino models, with all parameters in the superpotentials that are 
real. We will limit ourselves to the wall solutions where all fields take real values, 
so we do not have to travel in the complex plane, and can apply a rich physical 
intuition stemming from the fact that the BPS equations in this case are those of 
high- viscosity fluid (the so called creek equations) 0]. We will see that even in this 
simplest case the situation is quite non-trivial. Whenever the low-energy theory has 
a glued potential, we can count the number of distinct walls but, generically, can 
say nothing about their BPS nature and/or tension. 

Let us consider for simplicity two chiral superfields, $ and X, and a superpo- 
tential shown on Fig. 3. (More exactly. Fig. 3 displays W as a function of 0, x ^or 
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Figure 3: A superpotential with two mountain ridges and a canyon and five vacuum 
states. 



real values of Shown are two "mountain ridges", the left ridge and the right 

one, separated by a "canyon" . The heavy field is $, the light one is X. The vacua 
of the theory correspond to the points where 

These points are denoted by A, B, C, D, E on Fig. 4. The points A, B, C lie on the 
left ridge, the point D on the right ridge while the point E belongs to the bottom 
of the canyon. 

The low-energy reduction is obtained by eliminating the field $ by virtue of the 
equation 

'-^-0 (60) 

Substituting back in the superpotential W we get an effective low-energy 

superpotential 

Wefr(X) = (61) 
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Figure 4: The projection of the superpotential of Fig. 3 onto the Wx plane. Shown 
are the shapes of the mountain ridges and the canyon bottom. The points of extrema 
in W are denoted by A,B,C,D,E. 
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Equation ( |60D determines the positions of the mountain ridges and the bottom of 
the canyon, while Eq. (^) projects them onto the Wx plane. One can visualize 
Weff as shadows of the ridges left by a parallel horizontal beam of light in the (j) 
direction. 

First of all, let us prove that the vacua of the full theory (extrema of W) lie on 
the mountain ridges and/or bottoms of the canyons (extrema of Wcs )• Indeed, 



cff 



dX 



9W($,X) ^ 9W($,X) 9$, 



dX a$ dX 



= 0, 



if Eq. ( |59D is satisfied. Thus, the extrema of Wefr and W coincide. 

Let us assume now, for a short while, that the right mountain ridge and the 
canyon do not exist, and the profile of W has only the left mountain ridge, with 
three vacua. A, B and C . In this case the solution of Eq. (^) is unique. The full 
theory has two BPS walls: AB and CB. The potential of the low-energy effective 
theory >VefT(X) is a smooth function, and the low-energy theory has two BPS walls 
too. Their existence can be seen from the creek equation in the low-energy theory 
per se, 

d.X = ^. (62) 



The wall profile X{z) we find from Eq. (|6^) is a little bit different from what one 



would get by solving the creek equations in the full theory. 

The difference vanishes in the limit when the field $ becomes infinitely heavy; it 
dies off as positive powers of At the same time , the wall tension found in 

the effective theory exactly coincides with that one would find in the full theory. No 
corrections can appear. This is a remarkable feature of BPS supersymmetric 
walls. The tension of such a wall is exactly determined by the central charge 0, 
H, 0] which reduces, in turn, to |W($o,-^o) - >V($i,Xi)| in the full theory, and to 
|>V(Xo)efr — W(Xi)efr| in the effective low-energy theory. (Here {$0,15-^0,1} are the 
points of extrema.) The two expressions above coincide identically. 

Summarizing, if for all values of the light fields the solution for the heavy fields, 
to be integrated out, is unique, the existence of a BPS wall in the effective theory 
entails the existence of such a wall in the full theory, and vice versa. Moreover, 
if one calculates the wall tension in the effective theory, one gets the exact answer 
valid in the full theory, with no corrections. 

Let us return to the superpotential depicted on Fig. 3, with two mountain ridges 
and one canyon. In addition to the AB and CB walls, the full theory has a BE 
wall which is BPS saturated, and three continuous sets of BPS walls AE, CE and 
DE. Each set includes an infinite amount of degenerate walls (by degenerate we 



26 



mean that the tensions of all walls inside each set are the same). The phenomenon 



of continuously degenerate supersymmetric walls was first observed in Ref. ^5 
Besides these BPS walls, the full theory may have AD, BD and CD walls that are 
not BPS saturated. Depending on the values of parameters in the superpotential it 
may be expedient for some or all non-BPS walls to decay into a pair of BPS walls. 

What can be said about the effective theory? The corresponding low-energy 
effective potential will be glued out of five pieces, as indicated in Fig. 4. In each of 
five domains, the branch of the effective superpotential corresponding to the lowest 
energy is shown by rectangles. We get a typical five-sector structure of the scalar 
potential in the effective theory. /^From this structure we conclude, with certainty, 
that the theory under consideration has five degenerate vacuum states, so that each 
pair can be connected by a domain wall. Inspecting the low-energy theory, without 
additional information on the full theory, we can say nothing, generally speaking, 
as to the nature of these walls. The reason is obvious: we have no idea where each 
sector comes from, whether or not two extrema in question belong to one and the 
same mountain ridge (canyon). If they do, no restructuring of the heavy field (f) 
vacuum occurs inside the wall, and we return back to the situation with the unique 
solution of Eq. (^) discussed above. 

If extrema from different sectors actually do belong to distinct mountain ridges 
or canyons | (Eq. (|60|) has more than one solution) there is no unambiguous way 
to decide BPS versus non-BPS from inspection of the low-energy theory alone. We 
need additional information regarding what happens with the heavy fields inside the 
wall. In any case, a part of the wall tension associated with the light fields will not 
saturate e. 



4 The Kaluza-Klein Domain Wall 

As was discussed above, the tension of the walls interpolating between the neigh- 
boring vacua in supersymmetric gluodynamics is expected to scale as N rather than 
N"^. This gives rise to a natural identification of these walls with the D-brane so- 
lutions found in Ref. ||^. Here we will show, that this phenomenon, "abnormal" 
N dependence, is actually more general, and shows up in other wall configurations 
related to Witten's analysis. In fact, the low-energy limit of the theory considered 
in Ref. is a five-dimensional Kaluza-Klein (KK) theory with a five- dimensional 
SU{N) gauge field Am, M = 1,..,5 and charged matter. After compactification 
there are two types of gauge fields - our original SU{N) gauge field A^ and a new 
f/(l) gauge field -B^ coming from the G^s components of the metric tensor as well 
as a scalar made from the fifth component A^ of the SU{N) gauge field. We are 
going to demonstrate that if this theory is modified, so that the supersymmetry is 
broken explicitly, there is a new type of domain wall due to the field A^. By ana- 

^Generally speaking, distinct sectors can belong to one and the same branch, see points A and 
C in Fig. 4. 
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lyzing this low-energy theory per se, with no reference to D-branes, we demonstrate 
that the wall tension scales as A^, in parallel with the brane-based derivation of Ref. 
0. Moreover, these walls carry an induced fractional charge [0]. Conceptually the 
situation reminds that with the axion wall discussed in Sect. 2.2. 

The theory to be considered is gravity plus the gauge field Am in five dimensions. 
To warm up we start from the Abelian case, i.e. the U{1) gauge group. At this stage 
we also omit the superpartners from the discussion. The action is 



S = J d^xV^ 
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(64) 



where ^5 and 65 are the five-dimensional gravitational and gauge coupling constants, 
Gmn is the metric, R stands for the curvature, F^'^^ is the gauge field strength 
tensor, all capital Latin letters run from 1 to 5, say, M, N = 1, 5, while the Greek 
letters /i, i/ , ... = 1, 4. It is assumed that one of the five dimensions forms a circle, 
so that we deal with x KK model. The mater sector consists of charged 
scalars $, the simplest possible choice. 

After x decomposition of the metric 

C^MN - ^ ^ j^^j^, j,G - ^ j (65) 

we get four-dimensional gravity, two U{1) gauge fields and S^, plus a scalar 
(we put dilaton G^^ = 1). 

For any manifold K with a nontrivial iTilK) the KK theory contains special 
Wilson line operators 

= P exp (^i j> Adx^ , (66) 

where 7 is a closed noncontractible contour on K. In our case K = x 5*^ and 
7ri(S'^) = Z. For the U{1) gauge field 

U = e'^, (j) G [0, 27r) . 

The phase represents the constant component of the gauge field A^ = [0,1/ R], 
where R is the compactification radius. Other values of A^ can be gauged into the 
interval [0, 1/-R) while ^5 = is gauge equivalent to A5 = 0. The corresponding 
gauge function is e = x^/ R. This means that the charged (with respect to field A) 
fields are related by a gauge transformation, 

^A,=i/Rix, X5) = e'"^/^^A,=o(a;, X5) (67) 

where x denotes four noncompact coordinates. Using the standard KK decomposi- 
tion 

00 

vl>(a:,X5)= E e*""^/^vl>„(x) , (68) 
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where \l/n(a^) describes different four-dimensional fields with masses = v? jB? 
and charges g = n with respect to the KK gauge field = G^s, we see that the 
gauge transformation ( p^ ) shifts n ^ n + 1: 

^n{x)\A,=i/R = ^n+i(a;)Us=o (69) 

Changing adiabatically A<^ from zero to all levels in the particle spectrum 
are shifted by one. For example, a massless neutral (with respect to B^) particle will 
be transmuted into a heavy {m = 1/R) charged {q = 1) particle. Thus, if there is a 
domain wall ^5(0;) interpolating between y45(— 00) = and ^5(00) = 1/R, and one 
scatters a massless neutral particle with n = on the domain wall it either reflects 
or becomes massive and charged. The total KK charge must be conserved in the 
process, since it is the gauge charge which generates a part of the general coordinate 
transformation, 6Gmn = Qm^n^x^x^) + dNeuix^x^) which for eM{,x,x^) = Sm,5^{x) 
is the f/(l) gauge transformation, = ^G^s = 9^e. 

The conservation of the KK charge is insured in a way very similar to the one 
discussed in Sect. 2.2. In the presence of the charged particles the domain wall 
itself acquires an induced KK charge q. The total charge of the domain wall plus 
the particle is conserved. The real process of the particle penetration through the 
domain wall looks as follows: the initial state is the KK neutral particle(g = n = 0) 
plus the charged domain wall, with the charge gj = +1/2. The final state (if the 
particle initially has momentum px > 1/-R) is the charged particle q = n = 1 plus 
the domain wall with the charge q/ = —1/2. The total charge +1/2 is conserved. 

These walls have a variety of interesting properties. Thus, for example, if the 
theory contains fermions, their charge may be half integer. In this case the charges 
of the wall and the fermion are exchanged in the scattering process and there is no 
threshold energy for this process. Moreover the fermion and the wall can form a 



neutral bound state. For a detailed discussion see Ref. ||T8[| . 

Under what circumstances is this wall stable? The five-dimensional Maxwell 
term {l/AeDFl^^ gives rise to the four-dimensional kinetic term (1/46^)9^^59^745 
where e is the four-dimensional gauge coupling. It is quite evident that to get a 
stable domain wall an effective potential (which does not exist at the classical level) 
must be generated. Such an effective potential ^7(^5) is indeed induced at the 
one-loop level. The mass operator for the scalar field modes is 

rh = -185 + A5; mn = — + A5. (70) 

K 

The mass spectrum depends on the value of with periodicity 1/R. Due to this 
dependence there is an effective potential for the field A^ which is given, at one loop, 
by (see |^ for details) 



V{A,) = ^^^^|4^g^2 / 1^(1 - 2e-* cos 2rrRA, + e'^*) . (71) 



^ We refer to the charge q as the KK charge, as opposed to the conventional charge defined 
with respect to A^. 
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This potential is periodic in A5, with the period 1/R. It has minima sX A^ = n/R. 
Using this effective potential it is not difficult to see that the wall width L is of the 
order of R/e and the energy density (wall tension) scales as 

L{l/e'){l/RLf = {l/e){l/R% 

Observe that the wall tension scales as 1/e rather than 1/e^. This is essentially 
the same difference as between the and A^^ scaling laws in SUSY gluodynamics. 
Indeed, is an effective coupling constant in the string description of the gauge 
theory. 

The same N dependence arises also in the non-Abelian case. If there are no fields 
transforming in the fundamental representation, the gauge group G is actually not 
SU{N) but SU{N)/Zjq and 7ri(G') = Z^, so one may have a domain wall. 
The domain walls will be somewhat more complicated: instead of one phase field 
there are — 1 independent phases (=rank of the group SU{N)). 

In the supersymmetric theory the situation is essentially different. The fermions 
and the bosons give contributions of the opposite signs to the effective potential, 
which, thus, cancel each other. This is so since supersymmetry should be main- 
tained at any value of and the vacuum energy must be zero. Therefore, in 
five-dimensional SUSY gauge theory (which is essentially a low-energy limit of Wit- 
ten's theory) there are no stable domain walls of the type just discussed. To get 
a non-zero one-loop effective potential supersymmetry must be broken explicitly. 
This can be achieved either by imposing different boundary conditions on bosons 
and fermions, or by adding non-SUSY mass terms by hand, or in any other way 
- in any case there will be an induced effective potential. In the particular case of 
X U{1) relevant details can be found in Ref. [^, and the result for the wall 
tension in the large N limit is 



At large the proper coupling constant is not e but g = eyN, and the domain 
wall tension has a typical D-brane behavior - linear in both A^ and !/(?, 

e^{N/g){l/R'). 

Thus we see that in the theory of the type considered in (with explicitly 
broken SUSY) a new type of the domain wall, with an abnormal A^-dependence, 
emerges. The I^-brane interpretation of these walls should also be possible. 



"'^^Let us note that these walls are nothing but a 4 + 1-dimensional generalization of Z^r domain 
walls in gluodynamics at finite temperature p^ . The only difference is that for finite-T theory 
we must have anti-periodic boundary conditions for fermions, while in this case there is no such 
requirement. 
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5 Comment on Discrete Anomaly Matching 



The issue of the discrete anomaly matching in SUSY gluodynamics was raised in 



Ref . [|10l , with the conclusion that the chirally symmetric vacuum suggested in Ref . 
IP does not satisfy the matching condition, at least in its naive realization. Although 
this topic is peripheral to the main subject of the present work, we would like to 
dwell briefly on this issue in view of general important consequences which might 
follow. 

The construction which goes under the name of discrete anomaly matching was 



suggested in Refs. [^, |2^, 0, in a model-building context where it was quite 
informative and useful. Below we will show that in supersymmetric gluodynamics 
and similar settings no new physical results can be obtained from the procedure, 
besides those Ward identities which are already obtained by a different method. 
These Ward identities neither support nor rule out the chirally symmetric Kovner- 
Shifman vacuum. 

Let us introduce the construction in an explicit form, first in a somewhat simpli- 
fied setting of a non-supersymmetric Yang-Mills theory, with the intention to revisit 
SUSY gluodynamics later on. Assume we have a SU{N) Yang- Mills theory with 
one quark field belonging to the adjoint representation of the gauge group. Note 
that, unlike SUSY gluodynamics, \1/ is the Dirac field, so it can be coupled to an 
external "electromagnetic" field vectorially, with a very small coupling constant 
e, 

AC = e^' V^"^/. • (73) 

The field is auxiliary and is needed only for the purpose of constructing the 't 
Hooft AVV triangles [0. The quark is assumed to be massless. 

This theory, classically, has two conserved currents, the vector current ^'j'^'^ 
and the axial one, ^'jf^'-f^'^ . (There are two other conserved currents [0, but this is 
another story.) The vector current is useless from the point of view of the 't Hooft 
matching [^] and we will forget about it, focusing on the axial U{1) symmetry 
associated with the axial current. This symmetry is internally anomalous, 

d'^ (^7^'^) = £-,G%G;^ , (74) 

so that the only remnant of this "non-symmetry" is a discrete Z^n. The fact of 
survival of the discrete Z^n symmetry is most readily seen from the instanton- 
induced 't Hooft interaction ||33| which in the case at hand includes 4N fermion 



lines. The factor AN is related to the coefficient in the right-hand side of Eq. 
([T^), which in turn represents T{G), (one half) of the Dynkin index for the adjoint 
representation. If instead of the adjoint quark we were dealing with the quark in a 
representation R, then T{G) — > T{R). 

The survival of a discrete unbroken subgroup of the axial U{1) is similar to what 
we have in SUSY gluodynamics. 



31 



To construct the 't Hooft triangles that must be matched at the fundamental 
and constituent levels we have to have a continuous axial symmetry, rather than a 
discrete one. The basic idea of the discrete matching |]2S|, |2^, is embedding the 



theory under consideration into a larger one, where we do have a continuous axial 
symmetry, which is later spontaneously broken down to Z^n. The spontaneous 
breaking should happen in such a way that at low energies we recover the original 
theory plus possible extra decoupled degrees of freedom. Such an embedding is easy 
to achieve by exploiting the phantom axion construction [^. Let us add to our 



Yang-Mills theory a quark field Q in the fundamental representation of the gauge 
group. This quark is coupled to a (color-singlet) complex scalar field 

AU = hQRQL^ + }i.c. + Vm), (75) 

where h is a coupling constant, and y(|$|) is a self-interaction potential which will 
eventually ensure the development of a large vacuum expectation value of the field 

{^) = V ^ oo . 

When V becomes very large, the quark Q and the modulus of the filed $ disappear 
from the spectrum, leaving only Arg$, the axion field, as a remnant. The quark Q 
and the modulus of the filed $ are auxiliary elements of the construction. With the 
extra fields introduced in this way we have an additional axial U{1) symmetry 

Ql-^Qlc'", Qb -^Qrc-'^, $^$e-2^^ (76) 

This symmetry is internally anomalous too. Out of two internally anomalous U{1) 
symmetries we can readily pick up an anomaly-free combination. 

The corresponding conserved axial current has the form 

= i|r^/^^5^ _ 2NQYl^Q + scalar term . (78) 

The vacuum expectation value of the field $ spontaneously breaks the U{1) sym- 
metry of Eq. ([TTD, but since the corresponding charge of the $ field is 4iV, there is 
a survivor, a Z4n- subgroup. 

So, we have a theory which, at the fundamental level, has an internally anomaly- 
free axial U{1) containing a discrete subgroup. At the scale below v only the discrete 
subgroup survives. Let us examine now implications of the 't Hooft anomaly match- 
ing. The triangle to be considered is AVV where the A vertex is due to the axial 
current (|78|) , while the vector vertices are due to Eq. (|73|) . Note that the auxiliary 
quark Q has no coupling to A^. 

The AVV triangle appearing at the fundamental level is 

d'J, = ^-^F,^F,^ (79) 
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where F^^ is the photon field strength tensor built of the auxiliary field A^. If the 
photons are on mass shell, Eq. (^) implies ^ the existence of a pole coupled 



to J/,, 

(-^m) = ^ g^2 ^Mi'^Mi' • (80) 

The coefficient in front of g^/g^ has to be matched by the contribution of physical 
massless particles. Some of them may or may not occur dynamically, as composite 
mesons or baryons built from ^^'s in the original Yang-Mills theory under considera- 
tion. More important is the occurrence of the massless axion field, which is coupled 
to the current and, thus, participates in the matching with necessity. This is 
a distinctive feature of the discrete matching, as opposed to the 't Hooft match- 
ing, where such field, totally foreign to the original Yang-Mills theory per se, does 
not emerge. It is to be stressed that, as opposed to the Peccei-Quinn construction 



35|, in the present setup the axion is necessarily massless and can not acquire mass 
through nonperturbative effects. This is so since it is a Goldstone boson appearing 
due to spontaneous breaking of a continuous global symmetry. 

The axion field a is not coupled to the field because the auxiliary quark Q 
does not have this coupling. It is coupled, however, to the gluon field, through the 
standard vertex 

a, , . (81) 

Since its coupling to the current is 

jaxion _ _4]Vt;2(5^a) (82) 

we conclude that at low energies, in the effective low-energy theory, ( J^) = {J^)^^°^ 
plus a possible pole term in (J^) due to the contribution of massless composites 
built of should they exist. Here 

W'"'' = ^ (-^) (0|G'«.G«J77) • (83) 

The momentum q in Eq. (^) is the momentum flowing in the GG vertex (the total 
momentum of the photon pair). It is assumed that q —>■ 0. 
The matching of Eqs. (|0]) and (|83|) tells us that 

+ (^) {0\G%G';,M (84) 

= possible contributions due to massless composites. 
If there are none, then the expression on the left-hand side vanishes. 
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Figure 5: The two-photon matrix element of NGG. The photons are assumed to 
be on mass shell. 

Recall that our initial task was getting information on the emergence or none- 
mergence of the massless composites. The entire construction with embedding 
the discrete remnant of the anomalous axial U{1) was designed for that purpose. 
We are neither closer nor further now from this goal. Indeed, one can discard this 
construction altogether, and just consider the internally anomalous current (0). 
Then, combining both the external and internal anomaly, we would get 

{^ri'^) = ^,G%Gl^ + . (85) 

Sandwiching both sides of this formula between (0| and I77) in the limit g — ^ we 
immediately reproduce Eq. (P^). The only interesting dynamical question is whether 
the left-hand side of Eq. (|8^) vanishes or not. At first sight the dependence of 
two terms in this equation is different, so one is tempted to say that they cannot 
cancel each other. A closer look shows, however, that the discrepancy is superficial. 
Indeed, a typical graph for the second term is depicted in Fig. 5. The gluons are 
converted into photons through the ^ loop. It is not difficult to count that the 
matrix element shown in Fig. 5 scales as A^^, i.e. in the same way as the first term 
in Eq. (H). 

If T{R) were a free parameter than one could establish the non- vanishing of Eq. 
(p^) since the second term is proportional to (T(i?))^ while the first term to dimi?. 
The choice of T{R) is not free, however, since, on the one hand, to have a discrete 
unbroken subgroup Zj^ we must work with the quarks in representation higher than 
fundamental, but on the other hand the representation can not be too high, since 
otherwise we loose asymptotic freedom. For this reason T{R) cannot scale faster 
than A^. These two requirements are contradictory unless T{R) oc A^. 



Thus, Eq. (|8^) may or may not vanish, depending on whether the two terms 
cancel each other. As far as the A^ dependence is concerned, they are perfectly fit to 
cancel. In the absence of massless composites they would he forced to cancel. This is 
nothing but the NSVZ low-energy theorem for the two-photon coupling to GG 

Instead of the auxiliary photon we could have considered the coupling to 
gravitons, i.e. the U{1) current in the gravitational background. Then the issue 
would reduce to a formula connecting (A^^ — 1)RR to a two-graviton matrix element 
of NGG. Again, the so called discrete matching would have nothing to say whether 
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or not the two terms combine to cancel each other (in the first case there are no 
massless \1/ composites while in the second they would have to be present to match 
the anomaly). 

Now, we can readily adapt our consideration to SUSY gluodynamics. Again, 
we could have built a "tower of discrete anomaly matching" by embedding the 
theory in a larger one where an internally nonanomalous axial U{1) current exists, 
with the subsequent spontaneous breaking of this U{1) down to Z2N, the actual 
symmetry of SUSY gluodynamics. As we have just demonstrated, this procedure is 
redundant. It would yield no more constrains or information compared to what one 
gets considering the axial current of gluinos from the very beginning. The gluino 
is described by the Majorana field, so we cannot couple it to the auxiliary photon 
(vector current). However, the anomaly in the gravitational background remains an 
open possibility. The relation to be analyzed is 

d, (Wa) = ^(0|G;:,G^J2 gravitons) + -^-^CRR, (86) 

where C is a known constant. The question to be answered is: are there massless 
composites built from gluons/gluinos? In the standard chirally asymmetric phase 
we expect none, while in the chirally symmetric vacuum of Kovner and Shifman a 
set of massless composites must exist. 

We see that, if at all, the massless composites of the Kovner-Shifman solution 
facilitate the anomaly matching. Indeed, in the chirally asymmetric vacua the exact 
cancellation of two terms in Eq. (|86D must take place, while in the chirally symmetric 
one this cancellation can be partial. The missing part will then be filled by the 
contribution of massless composites. Regardless, the N dependence of both terms 
in Eq. (|86D is the same, and no constraints on the chirally symmetric solution |^ 
followQ. 

To summarize, the addition of extra matter which promotes the discrete Z2N 
symmetry into the continuous one necessarily leads to the appearance of a massless 
axion. This axion is indeed practically decoupled from the dynamics of the rest of 
the low energy sector when the symmetry breaking scale is large. However, at 
the same time it couples strongly to the conserved global current by virtue of Eq. 
(^). As a result, the contribution of the axion to the anomaly matching is finite 
and independent of the scale v. This contribution is the first term on the right hand 



side of Eqs. (|85|) and (|86|) . The discrete anomaly matching conditions, therefore. 



"'^"'^We note that the existence of the chirally symmetric phase was questioned recently on different 
grounds in [|6|. It was claimed that such a chirally symmetric phase would be necessarily super- 
conformally invariant and, therefore, have more symmetries than the Lagrangian of the original 
theory. Unfortunately this argument is not substantive. First, the fact of superconformal invari- 
ance of the chirally symmetric phase was not established in . It is perfectly conceivable that 
the correlators in this phase depend logarithmically on Aqcd- Second, even if the superconformal 
invariance is there, this is not forbidden by general principles of quantum field theory. For instance, 
in the realm of models of critical phenomena, the phenomenon of symmetry enhancement at the 
infrared fixed point is well known and not at all rare. 
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do not pose any restrictions on the spectrum of massless composite fermions but, 
rather, just determine the contribution of the axion which is not a physical quantity 
in the original theory. In this language the "modulo A^" matching of Ref. 



is the statement that the contribution of the axion to an anomaly triangle has to 
be an integer multiple Q of A^. This statement is true in a simple case when the 
fundamental fermions acquire mass only due to a Yukawa coupling to the scalar 
whose vacuum expectation value breaks the global U{1) symmetry down to the 
discrete subgroup in question. It is, however, not provable as a general result and 
there is no a priori reason to believe that it holds in strongly interacting theories 
with confining dynamics, like SUSY gluodynamics. 

For example, let us consider the very same toy model ([75|) , but this time, instead 
of the vector current ([73|), let us analyze the triangle with two vector currents Q'y^Q 
(and the same axial current as above). For convenience one can couple Q'J^Q to 
another auxiliary vector field Afj,, which is distinct from the field A introduced above. 
{A does not couple to ^.) The gauge strength tensor built from Afj, will be denoted 
by J-'fiiy- Now we can match the anomaly between the high energy scale M ^ v 
where the extra quarks Q are massless and the intermediate scale m ^ m ^ A). 
At this scale the quarks Q do not appear in the spectrum anymore. Note that the 
truly dynamical quarks \E' and their scale A are irrelevant in this problem. 

In this matching the contribution of the massless axion at the scale m must be 
equal to the contribution of the fermions Q at the scale M, since the only nontrivial 
dynamics that happens at v is the spontaneous symmetry breaking due to the Higgs 
field Since Q couple directly to $ their axial charges are necessarily integer 
multiples of 2N, Eq. ([77|). The contribution of Q to the anomalous triangle at M is 

1571-2 '"^^'^ 

the coefficient is multiple integer of 4N. Thus, it is indeed true that the contribution 
of the axion to this anomaly at the intermediate scale m is an integer multiple of 
AN. This is of course trivially so, since the axion coupling 

JLt t 

in conjunction with Eq. ( p2D automatically guarantee the required proportionality. 

We are interested, however, in matching the anomaly between m and a still lower 
scale fi < A. We would have the statement about "modulo N" matching here if we 
knew that the contribution of the axion changes by an integer multiple of 4N when 
crossing the scale A. This we can not know, however, because of the strong nontrivial 
interaction at A. The best we can do is to express this contribution in terms of a 



^^We use the phrase "integer multiple of TV" in a somewhat loose sense. This contribution 
depends on what particular anomaly triangle one considers and on some other details, e.g. the 
existence of massive Majorana fermions in the spectrum. These details are unimportant for us 
here. For a thorough discussion see Refs. [E9|, fSG] and [[lO|. 
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matrix element of GG, see Eq.(Q). The fact that the axion couples extremely weakly 
to the gauge field G does not help here, since it is a finite "renormalization" of this 
weak coupling due the strong interactions of the gauge field and \E' that determines 
the contribution of the axion to the anomaly. 

6 Discussion 

In this paper we have analyzed aspects of the supersymmetric walls in SUSY gluo- 
dynamics and in a more general context. We have argued that the linear scaling of 
the wall tension with does not contradict the picture where the wall is a classical 
soliton in an effective Lagrangian describing low-energy mesons and glueballs. We 
have also provided a simple qualitative explanation of how the confining string can 
end on a domain wall. 

A key part of our analysis is related to the issue what happens when the low- 
energy Lagrangian has a glued structure. The Veneziano-Yankielowicz description 
of SUSY gluodynamics belongs to this class. 

We have tried to answer the question whether the walls in SUSY gluodynamics 
are indeed BPS saturated. Our conclusion is that the knowledge of the effective 
Lagrangian by itself is not sufficient to answer this question if the wall in question 
crosses boundaries of distinct sectors. The effective Lagrangian has a cusp struc- 
ture which arises due to adiabatic integration of the heavy degrees of freedom. The 
adiabatic approximation breaks down when the wall trajectory crosses the cusp. 
We have considered a simple model which illustrates this feature in detail. In this 
toy model we were able to calculate the extra term, which gives the cusp contri- 
bution to the wall energy. Unfortunately in SUSY gluodynamics we are unaware 
of a well-defined procedure which could be used to obtain the missing term, since 
the (amended) Veneziano-Yankielowicz effective Lagrangian was not obtained by 
explicitly integrating out heavy fields, but, rather, from certain Ward identities of 
the theory. It is an interesting question whether these same Ward identities could 
also determine the cusp contribution. 

We then made some observations of the general nature pertinent to the theory of 
the domain walls in the effective low-energy theories. Supersymmetric walls possess 
some unique features. Namely, if we find the BPS wall tension (for a wall which 
does not cross the sector boundary) in the effective theory, the very same tension 
takes place in the full theory. There are no corrections inversely proportional to the 
masses of the heavy fields which were integrated out. 

Finally, we worked out the issue of the discrete anomaly matching in SUSY 
gluodynamics. This procedure, when appropriately implemented, is shown to impose 
no constraints on the existence of the chirally symmetric vacuum state ||^. 

We would like to make a remark on relation of our analysis to the calculations 
of P]. The analysis of P] is performed in the framework of the effective Taylor 
- Veneziano - Yankielowicz (TVY) effective Lagrangian, which in addition to the 
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"glueball" superfield S contains matter superfields corresponding to additional mat- 
ter fields in the SUSY QCD. Due to inclusion of these additional superfields the TVY 
Lagrangian does not have a cusp structure and the wall configurations considered 
in therefore do not cross any cusps. The numerical analysis of the simplest TVY 
Lagrangian carried in |g] showed that the BPS saturated wall, although present in 
the weak coupling regime - at small value of the Higgs mass - disappears for masses 
greater than some critical mass m*. Even more surprisingly at a slightly greater 
mass m^,* even the non BPS solutions that connect two adjacent vacua but do not 
pass through the chiral point S = disappear altogether. Since the VY effective 
Lagrangian is obtained from the TVY Lagrangian in the limit m — oo the authors 
of [H tentatively conclude that the only walls that exist in pure SUSY gluodynamics 
are the ones that pass through the chirally symmetric vacuum. If this is the case it 
is indeed very surprising, since the energy density of these walls in the large N limit 
is 0{N'^). We feel however that it is premature to draw such definite conclusions 
from the existing calculations for two main reasons. First, the calculations have only 
been performed at = 2, 3. It is not clear whether the critical mass remains finite 
at large N if the kinetic term in the TVY Lagrangian is taken to reproduce the 
correct large N scaling of meson masses. Secondly, even though the introduction 
of the matter fields eliminates the cusps in the effective potential we believe that 
in the SUSY gluodynamics limit this effect must be unphysical. The heavy degrees 
of freedom that should be responsible for smoothening of the effective potential in 
SUSY gluodynamics should be glueballs which are heavier than the order parameter 
field S but still much lighter than the matter fields, which become infinitely heavy 
and should completely decouple in this limit. To our mind it is therefore ques- 
tionable that the TVY effective Lagrangian reflects correctly the physics of heavy 
modes at large matter fields masses, where it ceases to be the effective Lagrangian in 
the Wilsonian sense. We think therefore that the existence or nonexistence of BPS 
saturated walls in SUSY gluodynamics remains an open question. We also would 



like to add that after this work has been finished we learned about the preprint ||37 
where the vacuum structure and domain walls in SUSY gluodynamics were studied 
using D-brane approach. 
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